Gate-induced ionization of single dopant atoms 
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Gate-induced wave function manipulation of a single dopant atom is a possible basis of atomic 
scale electronics. From this perspective, we analyzed the effect of a small nearby gate on a single 
dopant atom in a semiconductor up to field ionization. The dopant is modelled as a hydrogen-like 
impurity and the Schrodinger equation is solved by a variational method. We find that — depending 
on the separation of the dopant and the gate — the electron transfer is either gradual or abrupt, 
defining two distinctive regimes for the gate-induced ionization process. 

PACS numbers: 03.67.Lx, 85.30.De, 73.21.-b, 71. 55.-1 



The size regime where the discreteness of doping must 
be taken into account is brought within experimental 
reach by today's semiconductor lithography techniques. 
In this regime, single dopant atoms have been demon- 
strated to dominate the behavior of downscaled versions 
of conventional devices 0. On the other hand, the 
promising opportunity is offered to study the physics 
of semiconductors on their ultimate length scale by ad- 
dressing separate dopants. Putting a small gate close 
to a single impurity would, for example, allow for the 
manipulation of individual hydrogen-like wave functions. 
Furthermore, large electric fields (otherwise only achiev- 
able in astronomy) can be experimentally obtained in 
semiconductors due to the occurrence of large dielectric 
constants and small effective masses. Apart from the fun- 
damental importance, an ultimate application is found in 
a Si-based solid state quantum computer , in which 
the nuclear spins of single ^^P-dopants are envisioned as 
qubits. In this proposal, addressing a single qubit by 
NMR is achieved via the hyperfine interaction of the nu- 
clear spin and its valence electron, which can be tuned by 
modifying the electron wave function with a nearby gate. 
In a recent variation of this design 4], the ionization of 
single dopants by this gate is an essential ingredient. 

Our aim is to quantitatively investigate the effect of 
the electric field generated by a local gate on a single 
neutral dopant atom in a semiconductor, ultimately lead- 
ing to ionization. The response to small fields has been 
addressed before in the context of quantum computing 
0, Q . In this paper, the complete ionization process is 
discussed. Our approach incorporates the computation 
of time independent ground state wave functions of the 
system and, subsequently, the estimation of transition 
probabilities. We conclude that the separation of the 
dopant and the gate determines the nature of the ioniza- 
tion process. When the dopant resides close to the gate, 
the electron is gradually pulled away from the dopant 
when the gate voltage is increased, while for a larger sep- 
aration the dopant ionizes abruptly at a well-defined gate 
voltage. 

Addressing a single dopant requires a small local gate. 
When a dopant would be ionized by a large gate (e.g. 
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FIG. 1: The dashed line represents the calculated potential 
due to the gate at the symmetry axis of the device for ta = 
2 a.u. and a gate voltage of 2 a.u. The solid line includes 
the dopant potential for d = 10 a.u. (Note that e.g. in silicon 
1 a.u. ~ 3 nm for lengths and 1 a.u. ~ 90 mV for voltages.) 
The inset shows a schematic of the device layout, indicating 
the important parameters. 



an infinite strip i3j), the electron would be delocalized 
along the gate. This would be undesirable in applications 
where (spin-)phase coherence must be kept under control, 
such as a quantum computer. Therefore, we chose to 
model the gate as a circular disc, having the additional 
advantage that the complete system (dopant plus gate) 
is radially symmetric. The layout of our model system is 
schematically depicted in the inset of Fig. The disc- 
shaped metallic gate with radius ta is separated from the 
semiconductor bulk (relative dielectric constant Eg) by a 
barrier (relative dielectric constant Eb) of thickness dbar- 
A dopant is positioned at distance d from the barrier- 
semiconductor interface and centered with respect to the 
gate. 

At low temperatures, the semiconductor can be con- 
sidered as a dielectric, due to the absence of free charges. 
Charges at the barrier-semiconductor interfaces and in 
the barrier will be neglected. In our calculations, we as- 
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sume the barrier to be infinitely high and infinitely thin 
(<^bar = 0), which allows us to take advantage of the fact 
that the potential due to a charged metallic disc in a uni- 
form dielectric medium can be expressed in closed form 
(we will demonstrate the applicability of our results to 
a realistic layout). The total potential was obtained by 
adding a Coulomb potential well due to a positive unit 
charge in the semiconductor. A cross-section of the total 
electron potential for some typical parameters is shown in 
Fig. n Image charge effects at the semiconductor-barrier 
and the barrier-gate interfaces were neglected. 

In our calculations, the dopant atom is described 
within an effective mass approach: the contribution of 
the semiconductor bandstructure is accounted for by con- 
sidering it as an uniform dielectric medium and using 
an isotropic effective mass. Such a hydrogen-like model 
is known to provide a good first order description of a 
dopant atom (although it fails to accurately describe the 
energy levels Q and interactions 0). It is sufficient for 
our purpose and allows us to capture crucial phenomena 
and obtain estimates of important parameters. To keep 
our results general and transparent, physical quantities 
will be expressed in (effective) atomic units (a.u.) 
To simplify the conversion to conventional units, some 
values for silicon are given as an example in the caption 
of Fig. HI 

The time-independent Hamiltonian of the problem 
reads (in atomic units) 
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+ V^g(r,z), 



where V^{r,z) describes the potential landscape in the 
semiconductor due to the gate and (r, z) are cylinder 
coordinates as defined in Fig. ^ (inset). Approximate 
ground state wave functions are found by a variational 
method. As trial wave function we use a linear combina- 
tion of functions from a fixed and finite set S, where the 
weights are used as variational parameters. To this end, 
we choose S to contain functions of the form 



and 



(^(r, z) = cxp(— ar^) • z exp(— /3(z — dtf) (1) 



(^(r, z) = exp(— 7r^) • z exp(— (5z^), (2) 



where a, /3, 7, and 5 are constants that will be chosen 
later. The functions are cylinder-symmetric, motivated 
by the radial symmetry of the potential and the fact that 
the ground state is expected to be s-like. To allow for a 
full description of the ionization process, it is important 
that S includes both wave functions of the form , hav- 
ing large electron density at the dopant site, and of the 
form where the electron resides close to the gate. 

The functional form of Eq. is motivated by the 
fact that the (exponential) ground state wave function of 



hydrogen-like atoms can be quite well approximated as a 
linear combination of gaussians ,10.] , which are much eas- 
ier to work with numerically. To make sure that the wave 
functions vanish at the interface (2 = 0), it is multiplied 
by z. The (/?(r, z) are allowed to become aspherical due 
to the gate action by choosing different values for a and 
(3. Concerning the form of Eq. 10), we note that the po- 
tential well caused by the gate can in the radial direction 
be approximated by a parabola. Consequently, a ground 
state wave function similar to that of a linear harmonic 
oscillator is expected and therefore the r-dependent part 
of (p{r, z) is chosen as a gaussian. The ground state wave 
function of the triangular shaped well in the z-direction 
can be approximated as z • e~^^ (Ref. 0). Again, we will 
approximate the exponential by a linear combination of 
gaussians. 

In order to choose concrete values for the constants a, 
/3, 7, and 6, we note that for each positive integer N it 
is possible to find a set of N real numbers {Aij^^, such 
that a linear combination of exp(— A^r^) optimally ap- 
proximates the ground state wave function of hydrogen 
[lot . We wiU use the values given in Ref.0, which are, 
for example, {0.101, 0.321, 1.15, 5.06, 33.6} for iV = 5. In 
our calculation, we created functions of type by tak- 
ing values for a and (3 from such a set in all possible 
N"^ combinations. Functions of type were created 
by choosing values for 7 and 5 from the same set, after 

— 1/2 

multiplying all elements by the scaling constant to 
account for the size of the gate. Proceeding like this, S 
contains a total of 2N'^ functions. It was found that tak- 
ing N > 5 did hardly improve the accuracy. Therefore, 
iV = 5 was used in all presented results. 

Denoting the elements of 5 by the variational pro- 
cedure is now performed by forming the trial wave func- 
tion 

i^ir,^) = ^ Cn^n{r,z) 

as a linear combination of the and minimizing the 
functional 

m) 

with respect to the variational parameters c^. This min- 
imum is an upper bound to the ground state energy of 
Ti. This variational problem is equivalent to finding the 
smallest eigenvalue of the generalized matrix eigenvalue 
problem 



(H - £;M) -0 = 0, 



(3) 



where H is the hamiltonian matrix expanded on the ipn 
with elements Hij — ('0i|7i|i/>j) and M is the overlap 
matrix of the ipn defined as My = {tpi\i(;j). Further- 
more, c = (ci, C2, . . . , Cm) and the inner-product (-I-) is 



FIG. 2: (Color online) The radially integrated probability 
density of the electron wave function as a function of z for 
various gate voltages. The inset shows the ionization voltage 
versus gate radius ta. (In silicon: 1 a.u. ~ 3 nm and 90 mV, 
respectively.) 

(as usual) defined as 
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Note that M would be the unit matrix if S would be 
an orthonormal set with respect to (-I-). In that case, 
Eq. Q would reduce to an ordinary eigenvalue problem. 

The smallest eigenvalue Eq of Eq. Q is an upper 
bound to the ground state energy of the system. When S 
is chosen properly, Eq is a good approximation to the real 
ground state energy of H and the corresponding eigen- 
vector c defines a wave function that is a good approxi- 
mation of the real ground state wave function. 

Once this wave function is known for several values 
of the dopant depth d, gate voltage Vg and gate radius 
TA, we will use it to study the ionization process of the 
dopant. As an example, the radially integrated prob- 
ability density of the calculated electron wave function 
(i.e. \i{j{r, z)p27rr dr) is plotted versus z in Fig.EJ At 
zero gate voltage ^^e electron occupies the dopant 
site. For increasing gate voltage, the electron is gradu- 
ally pulled away from the donor site. Finally, for large 
enough gate voltage, it resides completely in the newly 
created potential well at the gate. 

An interesting physical quantity is the electron den- 
sity |?Ao(0, at the dopant site, as derived from the 
approximated ground state wave function ipQ(r,z). We 
will use it as an indication of the position of the electron: 
when the electron is pulled away from the dopant site, 
this number decreases. Moreover, it is of physical impor- 
tance because the hyperfine interaction is proportional to 
this number [l^ . 

The characteristics of the electron transfer from the 
dopant to the gate with increasing gate voltage depend on 
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FIG. 3: (Color online) The electron density at the dopant 
site as a function of gate voltage for various dopant depths, 
showing the process of ionization. All curves are normalized 
to their value at Vg = 0. The transition from a smooth to 
a step-like behavior is clearly visible at d ~ Sa^. The inset 
shows the corresponding data with a 2 a.u. thick oxide barrier 
present. The behavior is similar, but occurs at higher gate 
voltage. 

the distance d of the dopant under the gate. In Fig.O the 
electron density at the dopant site |'(/'o(0, (i)P (normalized 
to the value at zero gate voltage) is plotted as a function 
of gate voltage for several values of d. It can be seen that 
for small d the electron is transferred gradually from the 
dopant to the gate, while for larger d an abrupt electron 
jump occurs, defining an ionization voltage. This can 
be explained from the fact that for large d, a sufficiently 
large barrier separates the two potential wells. For small 
d, the two wells are so strongly coupled that they can be 
considered as a single well, the position of which is pulled 
towards the gate with increasing gate voltage. 

The calculations were repeated for several gate radii 
TA. The inset of Fig. [21 shows the ionization voltage for 
d = 15 versus ta. From the figure it is clear that the volt- 
age gets smaller for larger ta . The reason for this is that 
the transfer roughly takes place when the ground state 
energy of the gate-well drops below that of the dopant- 
well. When the gate-well is larger, the ground state en- 
ergy is closer to the bottom of the well and the transfer 
takes place at lower gate voltage. 

In a realistic device, the barrier between the gate and 
the semiconductor will have a finite thickness (in the 
most common material systems this will be at least 1 to 
2 a.u.). Usually, this barrier does not have the same di- 
electric constant as the semiconductor and hence it can 
modify the gate potential considerably. Moreover, the 
gate must be connected to the outside world by some 
kind of interconnect. Such an interconnect must be sep- 
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arated from the semiconductor by a much thicker barrier 
in order to sufficiently screen its potential There- 
fore, in a realistic device, the gate must be buried in a 
thick layer of barrier material. 

To allow for comparison with the idealized situation 
in which our calculations were carried out, several cal- 
culations were repeated with a realistic barrier present. 
To that end, we obtained the potential landscape due to 
the gate by solving the Poisson equation with a finite el- 
ement method (FEM) It was found that for typical 
realistic parameters (e.g. a Si02/Si-system with Sg = 12, 
Eb = 4, and dbar — 2 a.u.), the potential landscape in 
the semiconductor is qualitatively similar to the situa- 
tion where the gate is put directly on the semiconductor. 
As demonstrated in the inset of Fig. 01 the same phe- 
nomena are observed, but they occur at a higher voltage 
than in the absence of a barrier. The voltage drop over 
the barrier can roughly be accounted for by a linear scal- 
ing factor that depends on et, and dbar- Indeed, we find 
from the FEM-calculations that for the given parameters 
about 31% of the gate voltage drops in the semiconduc- 
tor. This number is similar to the observed ratio between 
the ionization voltages with and without a finite barrier 
thickness. This justifies the presentation of mainly re- 
sults obtained with an idealized barrier. 

As a final remark in our discussion of the barrier, we 
note that for any application or measurement of a single 
dopant device, it is crucial that there are no charge traps 
present near the dopant. Therefore it is highly desirable 
to have the barrier epitaxially grown on the semiconduc- 
tor. A promising candidate is a Sii-^^Gcx-layer as barrier 
on a Si substrate 0, although the maximum achievable 
barrier height in this system is only about 100 meV [l^ . 

The presented time-mdependent calculations are not 
sufficient to predict whether the dopant atom will indeed 
be ionized when the ground state wave function has a 
low electron density at the dopant site. In order to com- 
plete our analysis, an estimate of the tunnel probability 
is needed. This is obtained by comparison with the res- 
onance lifetime of a hydrogen atom in an electric field. 
The typical field strengths considered in the region be- 
tween the gate and the dopant site are very large (e.g. 
0.05-0.5 a.u. for ta = 2, d = 10, Vg = 2, see Fig^ . Us- 
ing a calculation of the Stark effect in hydrogen [l5j while 
taking the value of Ry* for silicon, it is found that the 
electron lifetime at the dopant site ranges roughly from 
0.1 ps to 1 ns. This can be interpreted as the time it 
takes for the dopant to be ionized when the gate voltage 
is switched on and justifies our interpretation of Fig. |21 
as the representation of an ionization process. 

Our general analysis can be readily applied, as we per- 
formed the calculations with parameters that are con- 
sistent with the quantum computer design mentioned. 
First, controlled tuning of the hyperfine interaction by 
the gate, which is required in Ref. z, is possible only when 
d is small enough: from Fig. |31 we estimate d < 6 a.u. 



Switching off the hyperfine interaction, as required in the 
'digital approach' can only be achieved for large sep- 
aration between dopant and gate {d > 10 a.u.). Hence, 
the dimensions of the device determine in which of both 
regimes operation takes place. Second, our analysis can 
be used to estimate the required gate voltage to tune the 
hyperfine interaction to a certain value (Fig. Third, 
it is found that the required voltage to fully ionize the 
dopant depends on ta, but it is nearly independent of d 
(Fig EI). 

In conclusion, we analyzed the wave function manipu- 
lation of a semiconductor dopant atom by a small elec- 
trostatic gate. We find that two regimes can be distin- 
guished for the ionization process of the dopant. For a 
dopant-gate separation smaller than ~ SaJ (e.g. ~ 24 nm 
for P in Si), the electron is gradually pulled out of the 
Coulomb potential of the dopant. When the dopant re- 
sides further away from the gate, the transfer takes place 
abruptly at a well-defined threshold field. Both regimes 
are accessible, since, e.g., epitaxial growth techniques al- 
low for sufficiently accurate positioning of the dopant un- 
der the gate. 
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